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ABSTRACT 

The Herschel Reference Survey is a guaranteed time Herschel key project and will be a benchmark 
study of dust in the nearby universe. The survey will complement a number of other Herschel key 
projects including large cosmological surveys that trace dust in the distant universe. We will use 
Herschel to produce images of a statistically-complete sample of 323 galaxies at 250, 350 and 500 /im. 
The sample is volume-limited, containing sources with distances between 15 and 25 Mpc and flux 
limits in the K-band to minimize the selection effects associated with dust and with young high-mass 
stars and to introduce a selection in stellar mass. The sample spans the whole range of morphological 
types (ellipticals to late-type spirals) and environments (from the field to the centre of the Virgo 
Cluster) and as such will be useful for other purposes than our own. We plan to use the survey to 
investigate (i) the dust content of galaxies as a function of Hubble type, stellar mass and environment, 
(ii) the connection between the dust content and composition and the other phases of the interstellar 
medium and (iii) the origin and evolution of dust in galaxies. In this paper, we describe the goals 
of the survey, the details of the sample and some of the auxiliary observing programs that we have 
started to collect complementary data. We also use the available multi-frequency data to carry out 
an analysis of the statistical properties of the sample. 

Subject headings: Galaxies: general, ISM; Infrared: galaxies; Submillimeter; Catalogs; Surveys 
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1. INTRODUCTION 

Understanding the processes that govern the formation 
and evolution of galaxies is one of the major challenges 
of modern astronomy. Ideally this work can be done by 
analyzing and comparing the physical, structural and 
kinematical properties of different objects at various 
epochs to model predictions. How primordial galaxies 
transform their huge gas reservoir into stars and become 
the objects that we now observe in the local universe is 
a key question. A wide coverage of the electromagnetic 
spectrum to probe the atomic (21 cm line) and molecular 
(generally through the 2.6 mm CO line) gas components, 
the different stellar populations (UV to near-IR spectral 
range), the dust (mid- and far-IR) and the electrons 
in magnetic fields (radio continuum) is necessary to 
address this important question. 

The importance of dust resides in the fact that it is 
formed from the aggregation of the metals produced 
in the latest phases of stellar evolution, and contains 
approximately half the metals in the interstellar medium 
(Whittct 1992). Injected into the interstellar medium 
(ISM) by stellar winds and supernovae explosions, dust 
acts as a catalyst in the process of transformation of 
the atomic to molecular hydrogen necessary to feed 
star formation (Hollenbach & Salpetcr 1971; Duley & 
Williams 1986). Dust also contributes to the shielding 
of the UV radiation field, preventing the dissociation 
of molecular clouds, and thus playing a major role in 
the energetic equilibrium of the ISM (Hollenbach & 
Tielens 1997). Dust contributes to the cooling and 
heating of the ISM in photodissociation regions through 
photoelectric heating. Furthermore, by absorbing the 
stellar light, dust modifies our view of the different 
stellar components (e.g. Buat & Xu 1996): because dust 
obscuration is so important in star-forming regions, the 
emission from dust is one of the most powerful tracers of 
the star formation activity in galaxies (Kennicutt 1998; 
Hirashita et al. 2003). Dust emission can therefore be 
used to study the relationship between the gas surface 
density and the star formation activity, generally known 
as the Schmidt law. 

The importance of dust in the study of the formation 
and evolution of galaxies became evident after the 
IRAS space mission which provided us with a whole 
sky coverage in four infrared bands. Despite its low 
sensitivity and poor angular resolution, IRAS detected 
tens thousands of extragalactic sources (Soifer et al. 
1987). The IRAS survey showed that the stellar light in 
some galaxies is so heavily obscured that the only way 
to determine their star formation activity is through 
the dust emission itself. Among these highly obscured 
galaxies, the Ultraluminous Infrared Galaxies (ULIRGs) 
are the most actively star-forming objects in the local 
universe. With their improved sensitivity, spectral and 
angular resolution, other space missions such as ISO 
and Spitzer have brought a better understanding of the 
chemical and physical properties of interstellar dust 
in a wide range of different galactic and extragalactic 
sources. IRAS, ISO, Spitzer and the recently launched 
AKARI satellite, however, are sensitive to dust emitting 
in the mid- (~ 5 /jm) to far- (< 200 /txm) infrared 



domain. The 5-70 /xm spectral range corresponds to the 
emission of the hot dust generally associated to star 
formation, while at longer wavelengths, up to 1 mm, 
the contribution of cold dust becomes more and more 
important (Draine et al. 2007). The emission of the ISM 
in the range 3 /xm < A < 15 /xm is generally dominated 
by Polycyclic Aromatic Hydrocarbons (PAHs). Between 
10 /im and < 70 /xm the dust emission is usually due to 
very small grains, while that at longer wavelengths (70 
A™ ^ $5 1000 M m ) i s probably produced by big grains 
of graphite and silicate in thermal equilibrium with the 
UV and optical photons of the interstellar radiation field 
(Desert et al. 1991; Dwek et al. 1997; Zubko et al. 2004; 
Draine & Li 2007). 

There is, however, strong empirical evidence to 
suggest that much of the dust in normal galaxies has 
been missed by previous space missions because it is too 
cold to radiate in the mid- and far-infrared. Devereux 
& Young (1990), for example, showed that when the 
dust masses of galaxies were estimated from IRAS data, 
the gas-to-dust ratios were ~10 times higher than the 
Galactic value, implying that 90% of the dust in galaxies 
is too cold to radiate significantly in the IRAS bands. 
ISO 200 /xm images of nearby galaxies revealed the 
presence of cold (^20 K) dust in the external parts of 
galaxies (Alton et al. 1998). Even Spitzer, with its 
longer wavelength coverage, is only sensitive to dust that 
is warmer than ~15 K (Bendo et al. 2003; Draine & Li 
2007; Draine et al. 2007). Simple estimates show that 
the temperature of a dust grain in thermal equilibrium 
in the average interstellar radiation field is ~ 20 K, 
which would produce a modified black body spectrum 
with a peak at ~ 200 /xm with a flux rapidly decreasing 
at longer wavelengths (Boselli et al. 2003; Dale et al. 
2005, 2007; Draine & Li 2007; Draine et al. 2007). It is 
important to remember, however, that most of the dust 
in a galaxy is likely to be at a much lower temperature 
than 20 K because it is self-shielded, and so unaffected 
from the interstellar radiation field. For example, in 
molecular clouds the temperature of the dust is usually 
well below the canonical 20 K value once the average 
dust extinction (Ay) is greater than one, unless it is 
close to a newly formed star (Ward-Thompson et al. 
2002). 

The submillimetre waveband (200-1000 fim) is crucial 
for detecting this missing cold dust component and thus 
for making accurate estimates of the total dust mass. 
Mass estimates made from far-infrared measurements are 
highly uncertain because of the strong dependence of flux 
on dust temperature which is difficult to determine using 
data on the Wien side of the black body peak. At long 
wavelengths (Rayleigh- Jeans domain) the emissivity de- 
pends only on the first power of the dust temperature, 
making it possible to use a submillimetre flux to make an 
accurate estimate of the mass of dust (Eales et al. 1989; 
Galliano et al. 2003; 2005). 

Submillimetre continuum observations of galaxies have 
been made previously, in particular with the ground- 
based submillimetre cameras, such as SCUBA on the 
James Clerk Maxwell Telescope. These have confirmed 
that galaxies do indeed contain a large amount of cold 
dust (Dunne & Eales 2001; Galliano et al. 2003). Until 
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now, the largest sub-millimetre survey has been the 
SCUBA Local Universe and Galaxy Survey (SLUGS), 
a survey of ~200 galaxies in two samples, one selected 
from the IRAS survey and one selected at optical 
wavelengths. The survey produced the first estimates of 
the submillimctre luminosity and dust-mass (the space 
density of galaxies as a function of dust mass) functions 
(Dunne et al. 2000; Vlahakis ct al. 2005), and also 
indicated that some ellipticals are surprisingly strong 
submillimetre sources. Its limitation, and the limitation 
of all ground-based submillimetre observations, is that 
of sensitivity: SLUGS struggled to detect galaxies not 
already detected by IRAS. 

The solution is to go into space and the Herschel space 
telescope is finally making this possible. In particular 
the relatively large field of view, high sensitivity, and 
coverage of a waveband (250-500 lira) in which galaxies 
are intrinsically bright, make SPIRE on Herschel the 
ideal instrument for a study of all extragalactic sources 
(Griffin et al. 2006, 2007; Wasket et al. 2007). SPIRE 
has the sensitivity to map large areas of the sky down 
to the confusion limit in quite modest integration 
times, providing submillimetre data for large samples of 
galaxies at different redshifs, and is thus well adapted for 
pointed observations for all kind of nearby extragalactic 
sources. 

The SPIRE consortium has defined a number of 
coordinated guaranteed time programs to get the best 
possible benefit of the unique Herschel facilities for the 
study of galaxy evolution. These include deep cosmo- 
logical surveys, pointed observations of local galaxies 
and surveys of representative regions in the Milky Way. 
To better characterize the dust properties in the local 
universe the SPIRE extragalactic group has selected a 
volume-limited sample of 323 galaxies to be observed 
in guaranteed time in the three SPIRE bands at 250, 
350 and 500 /im. The importance of the local universe 
resides in the fact that it represents the endpoint 
of galaxy evolution, providing important boundary 
conditions to models and simulations. Furthermore, at 
< 30 Mpc the angular resolution of SPIRE (a couple of 
kpc) allow us to resolve the different galaxy components 
such as nuclei, bulges, discs and spiral arms. Moreover, 
dwarf galaxies, by far the most common (yet still very 
poorly understood) galaxies in the universe, can only be 
observed locally. 

A volume-limited sample was chosen as a way to 
limit distance dependent selection biases. To limit any 
possible contamination due to the cosmic variance, the 
volume should be much larger than the typical dimen- 
sion of large-scale structures (~ 30 Mpc). At the same 
time it should be representative of the whole galaxy 
population inhabiting the local universe. A near-infrared 
K-band selection was adopted for two major reasons: 
(i) unlike optical surveys, which have strong selection 
effects due to the presence of dust, near-IR data are 
only marginally affected by dust extinction; (ii) whereas 
the optical flux is highly dependent on the number 
of relatively young stars and thus sensitive to recent 
episodes of star formation, the near-IR luminosity is a 
good measure of the overall stellar mass (e.g. Gavazzi 



ct al. 1996), which recent studies suggest as the most 
important parameter in characterizing the statistical 
and evolutionary properties of galaxies. Indeed galaxy 
properties that appear to be tightly correlated with the 
galaxy's stellar mass include the following: physical 
properties, such as the star formation activity, the gas 
content and the mctallicity (Boselli et al. 2001; 2002; 
Zaritsky et al. 1994; Gavazzi et al. 2004; Tremonti et al. 
2004); structural properties, such as the concentration 
index and the galaxy's light profile (Boselli ct al. 1997; 
Gavazzi et al. 2000; Scodeggio et al. 2002); the amount 
and distribution of dark matter, as indicated by the 
Tully-Fisher relation and the shape of the rotation 
curve for spirals (Tully & Fisher 1977; Catinella et 
al. 2006) and the fundamental plane for ellipticals 
(Dressier et al. 1987; Djorgovski & Davis 1987); the 
stellar population, shown through the colour-magnitude 
relations for early-type (Visvanathan & Sandage 1977, 
Bower et al. 1992) and late-type galaxies (Tully et 
al. 1982, Gavazzi et al. 1996) and the galaxy spectral 
energy distributions (Gavazzi et al. 2002, Kauffmann et 
al. 2003). These correlations all indicate a down-sizing 
effect (e.g. Gavazzi et al. 1996; Cowie et al. 1996; 
Heavens et al. 2004), in which galaxies with high 
stellar masses formed most of their stars at a much 
earlier cosmic epoch than those with low stellar masses. 
This underline the dominant role of mass, rather than 
morphology, in shaping galaxies. 

The final sample includes galaxies with a large range 
in luminosity (mass) and includes all Hubble types. Be- 
cause of its definition, the selected sample also includes 
galaxies belonging to different density regions such as 
the core of the Virgo cluster, groups and pairs and 
isolated objects. Given its completeness, the Herschel 
Reference Sample (HRS) will be a suitable reference for 
any statistical study. Combined with KINGFISH (the 
Herschel extension of SINGS, Kennicutt et al. 2003) and 
VNGS (see below), two samples optimized for the study 
of the different phases of the ISM in individual galaxies, 
these samples will provide the community with a unique 
dataset for studying the ISM properties of galaxies in 
the local universe. 

The paper is organized as follows. Section 2 describes 
the scientific goals of the survey. Section 3 gives a de- 
scription of the selection criteria used to define the sam- 
ple. Section 4 gives details of the SPIRE observations 
we will carry out and Section 5 gives an overview of 
the multi- frequency data that is available: for the sample:. 
Section 6 gives a brief description of the data processing 
and products. The multifrequency statistical properties 
of the sample will be described in Section 7. 

2. SCIENTIFIC OBJECTIVES 

The overall goal of this Herschel survey is to improve 
our knowledge of the cold dust properties of the most 
common extragalactic sources populating the local 
universe. Combined with multifrequency ancillary data 
covering the whole electromagnetic spectrum (see Sect. 
5), the new Herschel data will provide us with a unique 
dataset with which to: 

1) Trace, for the first time, the variation of 
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the properties of the cold dust component (dust 
mass and temperature, dust to gas ratio...) of 
normal galaxies along the Hubble sequence, 
and as a function of luminosity. Given the large 
dispersion in galaxy properties (Morton & Haynes 
1994), it is important the sample to be large enough 
to contain representatives of all galaxy types and 
include both early- and late-type galaxies spanning the 
largest possible range in luminosity. Where galaxies are 
resolved, we can analyze the distribution of the cold 
dust in the different galaxy components, e.g. the nuclei, 
bulges, spiral arms and diffuse discs. This will provide 
important constraints on dust formation and evolution 
in galaxies (Galliano et al. 2008). 

2) Study the role of dust in the physics of the 
ISM. As discussed in the introduction, dust plays a ma- 
jor role in the energetic equilibrium of the ISM. It acts 
as catalyst for the formation of the molecular hydrogen 
(Hollenbach & Salpeter 1971; Duley & Williams 1986) 
and shields the molecular gas component preventing 
dissociation from the diffuse interstellar radiation field 
(Hollenbach & Tielens 1997). To understand the nature 
of the ISM we thus need to know how the cold dust 
properties (temperature, composition, geometrical dis- 
tribution...) relate to other physical properties such as 
metallicity and interstellar radiation field (Boselli et al. 
2004; Galliano et al. 2008). By combining dust surface 
brightnesses with metallicity dependent gas to dust 
ratios and HI column densities, sub-mm measurements 
can be used to determine H2 column densities (Guclin 
et al. 1993; 1995; Neininger et al. 1996; Boselli et 
al. 2002). SPIRE data will therefore provide us with 
an independent measure of the molecular hydrogen 
component, and can be use for an accurate calibration 
of the still highly uncertain CO to H 2 conversion factor. 

3) Study the relationship between dust emission 
and star formation. Dust participates in the cooling 
of the gas through the shielding of the interstellar 
radiation field, in particular of the UV light, and thus 
plays a major role in the process of star formation 
(Draine 1978; Dwek 1986; Hollenbach & Tielens 1997). 
The energy absorbed by dust is then radiated in the 
infrared domain. For the same reason dust emission 
has often been used as a tracer of star formation. We 
still do not know, however, what is the relationship 
between the cold dust emission and star formation. The 
study of resolved galaxies will allow us to analyze the 
relationship between the infrared surface brightness and 
the dust temperature (Chanial et al. 2007) and to trace 
the connections between the star formation and the 
dust and gas properties at galactic scales, extending the 
recent results of Spitzer to all phases of the ISM (Gordon 
et al. 2004, Calzetti et al. 2005, 2007; Perez-Gonzalez 
et al. 2006; Kennicutt et al. 2007, Prescott et al. 2007; 
Thilker et al. 2007). 

4) Study the dust extinction properties in 
galaxies. Astronomers have tackled with the problem 
of dust opacity in galaxies for over 40 years (e.g. 
Holmberg 1958; Disney et al. 1989; Calzetti 2001), but 
have reached no consensus. The key problem is that 
optical catalogs contain huge selection effects because of 



the existence of dust. We will be able to address this 
issue in two ways. First, the submillimctre images will, 
for the first time, allow us to determine the distribution 
of the dust column density in a very large number of 
galaxies. Second, we will be able to determine the global 
dust opacity in each galaxy by using the energy balance 
between the absorbed stellar light and the dust emitted 
radiation (Buat & Xu 1996, Witt & Gordon 2000, Buat 
et al. 2002, Cortese et al. 2006; 2008a). This requires 
an accurate determination of the UV to near-IR (stellar 
light) and mid-IR to sub-mm (dust emission) spectral 
energy distribution (Boselli et al. 2003). By comparing 
the dust attenuation properties of different classes of 
objects this analysis will allow us to define standard 
recipes for correcting UV and optical data, a useful tool 
for the interpretation of all modern surveys. 

5) Determine whether there is an intergalactic 
dust cycle. Apart from the obvious possibility that 
dust is ejected from galaxies by the same methods that 
gas is ejected, such as through interactions with the sur- 
rounding environment (see point 9) and starburst-driven 
'superwinds' as indeed observed in M82 (Engelbracht 
et al. 2006), there is also the possibility that dust is 
ejected from galaxies by radiation pressure (Davies et 
al. 1998, Meiksin 2009, Oppenheimer & Dave 2008). 
The ejection of dust from galaxies might explain the 
huge reservoir of metals found in the intergalactic 
medium (Rayan- Weber et al. 2006). There is clear 
evidence for dust in superwinds (Alton et al. 1999) and 
there is tantalizing evidence from ISO observations for 
extended distributions of dust around galaxies (Alton 
et al. 1998; Davies et al. 1999). Observations so far 
have been limited by either sensitivity (SCUBA) or 
resolution (ISO and Spitzer). The Herschel Reference 
Survey will be able to determine whether dust ejection 
is important because (i) we will be observing several 
hundred galaxies, so even if these phenomena are rare 
our sample should contain some examples, and (ii) our 
observations will have the sensitivity to detect dust well 
outside the optical disc of each galaxy. 

6) Determine the amount of interstellar dust 
in ellipticals. Very little is known about dust in 
ellipticals. Despite the stereotype that ellipticals do not 
contain any dust, the structures seen in optical images 
imply that at least 50% of ellipticals contain some dust 
(Ferrarese et al. 2006). The amount of interstellar dust 
is too small, however, to be easily detected through 
its far infrared emission; IRAS, for example, detected 
only about 15% of ellipticals (Bregman et al. 1998). 
Although Spitzer has now detected many ellipticals, 
these studies have mostly been of sources known a 
priori to contain some dust (Kancda et al. 2007; Temi 
et al. 2007; Panuzzo et al. 2007) or molecular gas 
(Young et al. 2008). The SLUGS did contain a small 
statistically-complete sample of 11 ellipticals, although 
some of the six 850-^im detections may well have been of 
synchrotron rather than dust emission (Vlahakis et al. 
2005). Other ellipticals have been detected at 350 /im 
by Leeuw et al. (2008). The HRS contains 65 early-type 
(E, SO and SOa) galaxies and our observations will have 
sufficient sensitivity to detect dust masses down to 
~ 10 4 Mq. The Herschel Reference Survey will therefore 
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provide an unambiguous answer to the question of how 
much dust is in ellipticals. 

7) Determine the origin of dust in ellipticals. 

The three possible origins of the dust in ellipticals are 
that (i) it has been produced in the atmospheres of the 
old evolved stars that dominate the light of ellipticals 
today, (ii) its origin is external to the elliptical and 
is the result of a merger, (iii) it is the relic of dust 
formed during the active star-forming phase early in 
the history of the galaxy. The distribution of dust is 
an important clue to understand its true origin. The 
presence of the first mechanism is suggested by Spitzer 
mid-infrared spectra of Virgo ellipticals (Bressan et al. 
2006) which show the silicate emission produced by dust 
in circumstellar envelopes of evolved stars. The origin 
of this feature (predicted by Bressan et al. 1998) is 
supported by the finding that the mid-infrared emission 
has the same profile as optical light (Xilouris et al. 2004) 
in many early- type galaxies. It is unclear, however, 
whether this is the predominant mechanism producing 
the interstellar dust in these galaxies. For example, a 
submillimetre map of the nearby elliptical Centaurus 
A, has revealed a dusty disc, implying that the dust 
(and the galaxy itself) has been formed as the result of 
a merger (Leeuw et al. 2002). Temi et al. (2007) also 
found little correlation between the Spitzer 70 and 160 
/zm emission and optical starlight, which also suggests 
the dust has an external origin. 

We will use the Herschel Reference Survey results to 
investigate the origin of the dust by, first, investigating 
the detailed morphology of the dust, in particular how 
it compares with the stellar distribution, and, second, 
by investigating whether the mass of dust is correlated 
with other global properties of the galaxy such as stellar 
mass. It has often been argued that sputtering by the 
ubiquitous, hot X-ray emitting gas in early galaxies 
should destroy any dust formed more than 10 8 years in 
the past (Tielens et al. 1994). However, the fact that 
dust is seen visually in 50% of ellipticals (Ferrarese et 
al. 2006) and the recent discovery of small grains, which 
are preferentially destroyed by sputtering, in ellipticals 
(Kaneda et al. 2007) suggests that dust is protected 
in some way. We will investigate whether sputtering 
from the hot gas destroys the dust grains and will 
test the internal origin hypothesis by investigating any 
possible correlation between X-ray excess and dust mass. 

8) Measure the local luminosity and dust-mass 
distributions. These distributions are important 
benchmarks for the deep Herschel surveys. These have 
already been estimated as part of SLUGS (Dunne et al. 
2000; Vlahakis et al. 2005): the much greater sensitivity 
of the Herschel Reference Survey will push these limits 
down by a factor of ~ 50. 

9) Understand the role of the environment on 
the evolution of galaxies. The well-known phe- 
nomenon that spiral galaxies in clusters are Hi-deficient 
and have truncated HI discs (eg. Cayatte et al. 1990) 
demonstrates that the environment of a galaxy can have 
a strong effect on its interstellar medium. Indeed there 
is clear evidence for tidally-stripped dust in interacting 



galaxies (Thomas et al. 2002) and indications for 
ram-pressure stripped dust in cluster objects (Boselli 
& Gavazzi 2006). At the same time the presence of 
metals (Sarazin 1986) and possibly of dust (Stickel et 
al. 2002; Montier & Giard 2005) in the diffuse intra- 
cluster medium has been shown by X-ray and far-IR 
observations. The comparison of the submillimetre 
emission of cluster and isolated galaxies within the HRS 
will allow us to make a detailed study on the effects 
of the environment on the dust properties of galaxies, 
and thus understand whether the hot and dense cluster 
intergalactic medium can be polluted through the gas 
stripping process of late- type galaxies (Boselli & Gavazzi 
2006). 

10) Provide a multi-frequency reference sample 
suitable for any statistical study. Our aim is 

that the HRS will be the first large sample of galaxies 
with observations of all phases of the ISM, as well as 
measurements over the entire wavelength range of the 
spectral energy distribution (SED) for each galaxy. The 
sample will then serve for many purposes e.g. the UV 
to radio continuum SEDs could, for example, be used 
to determine the nature of unresolved sources or as 
templates for estimating photometric redshifts. 

3. THE SAMPLE 

The Herschel Reference Sample (HRS) has been 
selected according to the following criteria: 

1) Volume-limited: A volume limit was imposed 
to reduce distance uncertainties due to local peculiar 
motions and ensure the presence of low-luminosity, 
dwarf galaxies, not accessible at high redshift. By 
applying a lower distance (D) limit we exclude the very 
extended sources, the observing of which would be too 
time consuming We have selected all galaxies with 
an optical recessional velocity (vel, taken from NED) in 
between 1050 km s _1 and 1750 km s _1 that, for Hq = 70 
km s _1 Mpc -1 and in the absence of peculiar motions, 
corresponds to 15 < D < 25 Mpc E3- In the Virgo 
cluster region (12h< R.A.(2000) < 13h; 0°< dec < 18°), 
where peculiar motions are dominant, we have included 
all galaxies with vel < 3000 km s _1 and belonging to 
cluster A, the North (N) and East (E) clouds and the 
Southern extension (S) (17 Mpc) and Virgo B (23 Mpc), 
where the subgroup membership has been taken from 
Gavazzi et al. (1999a). W and M clouds objects, at a 
distance of 32 Mpc, have been excluded 

2) K band selection: Given the expected low dust 

23 A small sample of very nearby and extended galaxies will 
be observed in detail as part of another guaranteed time project: 
Physical Processes in the Interstellar Medium of Very Nearby 
Galaxies. 

24 Given the possible discrepancy between optical and HI reces- 
sional velocity measurements, heliocentric velocities given in Table 
1 can be outside this range. 

25 The sample was selected before distances were available on 
NED. According to these new NED estimates, the distance of the 
HRS galaxies outside the Virgo cluster are generally included in the 
15 < D < 25 Mpc range. In the Virgo region, however, our dis- 
tances are determined according to subgroup membership criteria, 
and are thus generally different than those given by NED 
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content of quiescent galaxies, whose emission would be 
hardly detectable within reasonable integration times, 
a more stringent limit has been adopted for early-types 
than for star forming galaxies. Among those galaxies in 
the required recessional velocity range, we thus selected 
all late-type spirals and irregulars (Sa-Sd-Im-BCD) with 
a 2MASS K band total magnitude Kstot < 12 and all 
ellipticals and lenticulars (E, SO, SOa) with Kstot < 8.7. 

3) High galactic latitude: To minimize galactic 
cirrus contamination, galaxies have been selected at 
high galactic latitude (b > + 55°) and in low galactic 
extinction regions (Ab < 0.2; Schlegel et al. 1998). 

The resulting sample is composed of 323 galaxies 
located in the sky region between 10hl7m < R.A.(2000) 
< 14h43m and -6° < dec < 60° (see Fig. [[]), of which 
65 are early-type (E, SO and SOa) and 258 are late-type 
(Sa-Sd-Im-BCD) (see Fig. ^. Figure H shows that 
the only galaxies which are clearly undersampled are 
blue compact galaxies and dwarf irregulars, the most 
numerous galaxies in the nearby universe. T hey will 
be the targets of another SPIRE key program ™ As 
selected, the sample spans a large range in environment 
since it includes the Virgo cluster, many galaxy groups 
and pairs as well as relatively isolated objects (Fig. [T]). 
Using the Virgo cluster membership criteria defined in 
Gavazzi et al. (1999a), the HRS includes 82 members 
of cluster A and B (Fig. [5]) The other galaxies are 
members of nearby clouds such as Leo, Ursa Major 
and Ursa Major Southern Spur, Crater, Coma I, Canes 
Venatici Spur and Canes Venatici - Camelopardalis 
and Virgo-Libra Clouds (Tully 1988). As defined, the 
sample is thus ideal for environmental studies (Cortese 
& Hughes 2009; Hughes & Cortese 2009). The whole 
sample is given in Table [T] with columns arranged as 
follow: 

Column 1: Herschel Reference Sample (HRS) name. 
This is the position of the galaxy in the sample list when 
sorted by increasing right ascension. 
Column 2: Zwicky name, from the Catalogue of Galax- 
ies and of Cluster of Galaxies (CGCG; Zwicky et al. 
1961-1968). 

Column 3: Virgo Cluster Catalogue (VCC) name, from 
Binggeli et al. (1985). 

Column 4: Uppsala General Catalog (UGC) name 
(Nilson 1973). 

Column 5: New General Catalogue (NGC) name. 
Column 6: Index Catalogue (IC) name. 
Columns 7 and 8: J2000 right ascension and declina- 
tions, taken from NED. 

Column 9: Morphological type, taken from NED, or 
from our own classification if not available. 
Column 10: Total K band magnitude (Kstot), from 
2MASS (Skrutskie et al. 2006). 

Column 11: Optical isophotal diameter (25 mag 
arcsec -2 ), from NED. 

Column 12: Heliocentric radial velocity (in km s _1 ), 
from HI data when available, otherwise from NED. 
Column 13: Cluster or cloud membership, from Gavazzi 



et al. (1999a) for Virgo and Tully (1988) or Nolthenius 
(1993) whenever available, or from our own estimate 
otherwise. 

Column 14: Pair/group membership, from Karachentsev 
et al. (1972) or from NED whenever available, or from 
our own estimate elsewhere. Pair/group membership 
has been assigned according to the following criteria: 
close pairs (CP) are those with a nearby companion at 
a distance less than 50 kpc and a difference in redshift 
< 600 km s -1 , as in Gavazzi et al. (1999b), while 
pairs (P) up to 150 kpc. The number indicates whether 
the galaxy belongs to a triplet (3), quadruplet (4) and 
quintuplet (5). Groups outside Virgo and its immediate 
surroundings have been determined by counting the 
number of bright galaxief(3 within 25 arcmin (which, 
at a median distance of 20 Mpc, corresponds to ~ 150 
kpc) and 600 km s . Pairs in the Virgo region are only 
those cataloged by Karachentsev et al. (1972). 
Column 15: Galactic extinction Ab (Schlegel et al. 
1998). 

Column 16: Alternative names. 



6 The ISM in Low Metallicity Environments: Bridging the Gap 
Between Local Universe and Primordial Galaxies 



27 For bright galaxies we intend those included in major catalogs 
such as NGC, UGC, IC or CGCG. 
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TABLE 1 
The Herschel Reference Sample. 



HRS 


CGCG 


VCCUGCNGCIC 


RA(2000) 
h m a 


dec 

O , ii 


type 


K S tot 
mag 


D(25) 


vcl 
km a" 1 


momb 


gror 


ipAB 
mag 


Other name 


1 


123035 


- 


- 


101739.66 


224835.9 


Pec 


11.59 


1 


00 


1175 


Leo Cl. 







13 




2 


124004 


5588 


- 


102057.13 


252153.4 


S? 


11.03 





52 


1291 


Loo Cl. 







10 




3 


94026 


5617 


3226 - 


102327.01 


195354.7 


E2:pcc;LINER;Sy3 


8.57 


3 


16 


1169 


Loo Cl. 


CP 





10 


KPG234A, ARP94 


4 


94028 


5620 


3227 - 


102330.58 


195154.2 


SAB(s)pec;Syl.5 


7.64 


5 


37 


1148 


Leo Cl. 


CP 





.10 


KPG234B, AR.P94 


5 


94052 


- 


610 


102628.37 


201341.5 


Sc 


9.94 


1 


86 


1170 


Leo CI. 







.09 




6 


154016 


5662 


3245A 


102701.16 


283821.9 


SB(s)b 


11.83 


3 


31 


1322 


Loo Cl. 


CP 





12 


NGC 3245A 


7 


154017 


5663 


3245 - 


102718.39 


283026.6 


SA(r)0:?;HII;LINER 


7.86 


3 


24 


1314 


Leo Cl. 


CP 





11 




S 


154020 


5685 


3254 - 


102919.92 


292929.2 


SA(s)bc;Sy2 


8.80 


5 


01 


1356 


Loo Cl. 







09 




9 


154026 


5731 


3277 - 


103255.45 


283042.2 


SA(r)ab;HII 


8.93 


1 


95 


1415 


Leo Cl. 







11 




10 


183028 


5738 


- 


103429.82 


351524.4 


S? 


11.31 





91 


1516 


Leo Cl. 







12 




11 


124038 


5742 


3287 - 


103447.31 


213854.0 


SB(s)d 


9.78 


2 


09 


1325 


Loo Cl. 







10 




12 


124041 


- 


- 


103542.07 


260733.7 


cl 


11.98 





59 


1392 


Loo Cl. 







10 




13 


183030 


5753 


3294 - 


103616.25 


371928.9 


SA(s)c 


8.38 


3 


55 


1573 


Leo Cl. 







08 




14 


124045 


5767 


3301 - 


103656.04 


215255.7 


(R')SB(rs)0/a 


8.52 


3 


55 


1341 


Leo Cl. 







10 




15 


65087 


5826 


3338 - 


104207.54 


134449.2 


SA(s)c 


8.13 


5 


89 


1300 


Leo Cl. 


3 





14 




16 


94116 


5842 


3346 - 


104338.91 


145218.7 


SB(rs)cd 


9.59 


2 


69 


1260 


Leo Cl. 







12 




17 


95019 


5887 


3370 - 


104704.05 


171625.3 


SA(s)c 


9.43 


3 


16 


1281 


Leo Cl. 







13 




18 


155015 


5906 


3380 - 


104812.17 


283606.5 


(R')SBa? 


9.92 


1 


70 


1604 


Loo Cl. 







11 




19 


184016 


5909 


3381 - 


104824.82 


344241.1 


SB pec 


10.32 


2 


04 


1630 


Leo Cl. 







09 




20 


184018 


5931 


3395 2613 


104950.11 


325858.3 


SAB(rs)cd pec: 


9.95 


2 


09 


1617 


Loo Cl. 


CP 





11 


KPG249A, ARP270 


21 


155028 


5958 


- 


105115.81 


275054.9 


Sbc 


11.56 


1 


45 


1182 


Leo Cl. 


3 





11 




22 


155029 


5959 


3414 - 


105116.23 


275830.0 


SO pec;LINER 


7.98 


3 


55 


1414 


Leo Cl. 


3 





11 




23 


184028 


5972 


3424 - 


105146.33 


325402.7 


SB(s)b:?;HII 


9.04 


2 


82 


1501 


Loo Cl. 


4 





10 




24 


184029 


5982 


3430 - 


105211.41 


325701.5 


SAB(rs)c 


8.90 


3 


98 


1585 


Leo Cl. 


4 





10 




25 


125013 


5995 


3437 - 


105235.75 


225602.9 


SAB(rs)c: 


8.88 


2 


51 


1277 


Leo Cl. 







08 




26 


184031 


5990 




105238.34 


342859.3 


Sab 


11.71 


1 


35 


1569 


Leo Cl. 







08 




27 


184034 


6001 


3442 - 


105308.11 


335437.3 


Sa? 


10.90 





62 


1734 


Loo Cl. 







.08 




28 


155035 


6023 


3451 - 


105420.86 


271422.9 


Sd 


10.23 


1 


70 


1332 


Leo Cl. 







09 




29 


95060 


6026 


3454 - 


105429.45 


172038.3 


SB(s)c? sp;HII 


10.67 


2 


.09 


1101 


Leo Cl. 


4 





15 


KPG257A 


30 


95062 


6028 


3455 - 


105431.07 


171704.7 


(R')SAB(rs)b 


10.39 


2 


38 


1105 


Leo Cl. 


4 





1-1 


KPG257B 


31 


267027 


6024 


3448 - 


105439.24 


541818.8 


10 


9.47 


5 


62 


1374 


Ursa Maj. S S 


CP 





05 


ARP205 


32 


95065 


6030 


3457 - 


105448.63 


173716.3 


S? 


9.64 





91 


1158 


Loo Cl. 


4 





13 




33 


95085 


6077 


3485 - 


110002.38 


145029.7 


SB(r)b: 


9.46 


2 


10 


1432 


Leo Cl. 







09 




34 


95097 


6116 


3501 - 


110247.32 


175922.2 


Scd 


9.41 


3 


89 


1130 


Leo Cl. 


P 





10 


KPG263A 


35 


267037 


6115 


3499 - 


110311.03 


561318.2 


10 


10.23 





81 


1522 


Ursa Maj. S S 







04 




36 


155049 


6118 


3504 - 


110311.21 


275821.0 


(R)SAB(s)ab;HII 


8.27 


2 


69 


1536 


Leo Cl. 


P 





12 




37 


155051 


6128 


3512 - 


110402.98 


280212.5 


SAB(r»)c 


9.65 


1 


62 


1373 


Loo Cl. 


P 





12 




38 


38129 


6167 


3526 - 


110656.63 


071026.1 


SAc pec sp 


10.69 


1 


91 


1419 


Leo CI. 







.14 




39 


66115 


6169 


- 


110703.35 


120336.2 


Sb: 


11.13 


1 


86 


1557 


Loo Cl. 







.07 




40 


67019 


6209 


3547 - 


110955.94 


104315.0 


Sb: 


10.44 


1 


91 


1584 


Leo Cl. 


P 





10 




41 


96011 


6267 


3592 - 


111427.25 


171536.5 


Sc? sp 


10.78 


1 


78 


1303 


Leo Cl. 







.07 




42 


96013 


6277 


3596 - 


111506.21 


144713.5 


SAB(rs)c 


8.70 


4 


06 


1193 


Leo Cl. 







10 




43 


96022 


6299 


3608 - 


111658.96 


180854.9 


E2;LINER: 


8.10 


3 


16 


1108 


Leo Cl. 


5 





09 


KPG278B 


44 


96026 


6320 


- 


111817.24 


185049.0 


S? 


10.99 





89 


1121 


Loo Cl. 


CP 





10 




45 


291054 


6330 


3619 - 


111921.60 


574527.8 


(R)SA(s)0+: 


8.58 


2 


69 


1544 


Ursa Major Cl. 


5 





08 




46 


96029 


6343 


3626 - 


112003.80 


182124.5 


(R)SA(rs)0+ 


8.16 


2 


69 


1494 


Leo Cl. 


CP 





09 




47 


156064 


6352 


3629 - 


112031.82 


265748.2 


SA(s)cd: 


10.50 


2 


29 


1507 


Leo Cl. 







08 




48 


268021 


6360 


3631 - 


112102.85 


531011.0 


SA(s)c 


7.99 


5 


01 


1155 


Ursa Major Cl. 







07 




49 


39130 


6368 


3640 - 


112106.85 


031405.4 


E3 


7.52 


3 


98 


1251 


Loo Cl. 


4 





19 




50 


96037 


6396 


3655 - 


112254.62 


163524.5 


SA(s)c:;HII 


8.83 


1 


55 


1500 


Leo Cl. 







.11 




51 


96038 


6405 


3659 - 


112345.49 


174906.8 


SB(s)m? 


10.28 


2 


09 


1299 


Leo Cl. 







08 




52 


268030 


6406 


3657 - 


112355.57 


525515.5 


SAB(rs)c pec 


10.29 


1 


45 


1204 


Ursa Major Cl. 







07 




53 


67071 


6420 


3666 - 


112426.07 


112032.0 


SA(rs)c: 


9.23 


-1 


37 


1060 


Leo Cl. 







14 




54 


96045 


6445 


3681 - 


112629.80 


165147.5 


SAB(r)bc;LINER 


9.79 


2 


25 


1244 


Loo Cl. 


4 





11 




55 


96047 


6453 


3684 - 


112711.18 


170149.0 


SA(rs)bc;HII 


9.28 


2 


.89 


1158 


Leo Cl. 


4 





11 




56 


291072 


6458 


3683 - 


112731.85 


565237.4 


SB(s)c?;HII 


8.67 


1 


86 


1708 


Ursa Major Cl. 


3 





07 




57 


96049 


6460 


3686 - 


112743.95 


171326.8 


SB(s)bc 


8.49 


3 


19 


1156 


Leo Cl. 


4 





10 




58 


96050 


6464 


3691 - 


112809.41 


165513.7 


SBb? 


10.51 


1 


35 


1067 


Leo Cl. 


4 





11 




59 


67084 


6474 


3692 - 


112824.01 


092427.5 


Sb;LINER;HII 


9.52 


3 


.16 


1717 


Leo Cl. 







1-1 




60 


268051 


6547 


3729 - 


113349.34 


530731.8 


SB(r)a pec 


8.73 


2 


82 


991 


Ursa Major Cl. 


P 





05 


KPG209B 


61 


292009 


6575 




113626.47 


581129.0 


Scd:;HII 


11.40 


1 


95 


1217 


Ursa Major Cl. 


4 





07 




62 


186012 


6577 


3755 - 


113633.37 


362437.2 


SAB(rs)c pec 


10.60 


3 


16 


1571 


Ursa Maj. S S 







.10 




63 


268063 


6579 


3756 - 


113648.02 


541736.8 


SAB(rs)bc 


8.78 


4 


17 


1289 


Ursa Major Cl. 







05 




64 


292017 


6629 


3795 - 


114006.84 


583647.2 


ScjHII 


10.64 


2 


14 


1213 


Ursa Major Cl. 


3 





06 




65 


292019 


6640 


3794 - 


114053.42 


561207.3 


SAB(s)d 


11.60 


2 


24 


1383 


Ursa Major Cl. 







06 




66 


186024 


6651 


3813 - 


114118.65 


363248.3 


SA(rs)b: 


8.86 


2 


24 


1468 


Ursa Maj. S S 







08 




67 


268076 


6706 




114414.83 


550205.9 


SB(s)m: 


11.28 


1 


91 


1436 


Ursa Major Cl. 







06 


NGC 3846A 


68 


186045 


- 




114625.96 


345109.2 


S? 


11.44 





32 


1412 


Ursa Maj. S S 







09 


MR.K 429 


69 


268088 


6787 


3898 - 


114915.37 


560503.7 


SA(s)ab;LINE;HII 


7.66 


4 


37 


1171 


Ursa Major Cl. 







09 




70 




- 


2969 115231.27 


-035220.1 


SB(r)bc?;HII 


11.15 


1 


23 


1617 


Crater Cl. 


3 





12 




71 


292042 


6860 


3945 - 


115313.73 


604032.0 


SB(rs)0 + ;LINER 


7.53 


5 


25 


1259 


Ursa Major Cl. 







12 




72 


- 


- 


3952 2972 


115340.63 


-035947.5 


IBm: sp;HII 


11.01 


1 


58 


1577 


Crater Cl. 


3 





11 




73 


269013 


6870 


3953 - 


115348.92 


521936.4 


SB(r)bc;HII/LINER 


7.05 


6 


92 


1050 


Ursa Major Cl. 


P 





.13 




74 


269019 


6918 


3982 - 


115628.10 


550730.6 


SAB(r)b:;HII;Sy2 


8.85 


2 


34 


1108 


Ursa Major Cl. 


4 





.06 




75 


269020 


6919 




115637.51 


553759.5 


Sdm: 


11.56 


1 


45 


1283 


Ursa Major Cl. 


4 





06 




76 


269022 


6923 




115649.43 


530937.3 


Im: 


11.32 


2 


00 


1069 


Ursa Major Cl. 


4 





12 




77 


13033 


6993 


4030 - 


120023.64 


-010600.0 


SA(s)bc;HII 


7.33 


4 


.17 


1458 


Crater Cl. 


P 





11 




78 


98019 


6995 


4032 - 


120032.82 


200426.0 


Im: 


10.45 


1 


86 


1269 


Coma I Cl. 







15 




79 


69024 


7001 


4019 755 


120110.39 


140616.2 


SBb? sp 


11.33 


2 


40 


1508 


Virgo Out. 







14 




80 


69027 


7002 


4037 - 


120123.67 


132403.7 


SB(rs)b: 


10.11 


2 


51 


932 


Virgo Out. 







12 




81 


13046 


7021 


4045 - 


120242.26 


015836.4 


SAB(r)a;HII 


8.75 


3 


00 


2011 


Virgo Out. 







10 




82 


98037 






120335.94 


160320.0 


Sab 


11.19 





60 


931 


Virgo Out. 







13 


KUG1201 + 163 


83 


41031 


7035 




120340.14 


023828.4 


SB(r)a:;HII 


11.82 


1 


.10 


1232 


Crater Cl. 







12 




84 


69036 


7048 


4067 - 


120411.55 


105115.8 


SA(s)b: 


9.90 


1 


20 


2424 


Virgo Out. 







11 




85 


243044 


7095 


4100 - 


120608.60 


493456.3 


(R')SA(rs)bc;HII 


8.03 


5 


37 


1072 


Ursa Major Cl. 







.10 




86 


41041 


7111 


4116 - 


120736.82 


024132.0 


SB(rs)dm 


10.27 


3 


80 


1309 


Virgo Out. 


P 





10 


KPG322A 


87 


69058 


7117 


4124 - 


120809.64 


102243.4 


SA(r)0 + 


8.49 


-1 


10 


1652 


Virgo Out. 







12 




88 


41042 


7116 


4123 - 


120811.11 


025241.8 


SB(r)c;Sbrst;HII 


8.79 


5 


00 


1326 


Virgo Out. 


P 





09 


KPG322B 


89 


69088 


66 7215 


4178 - 


121246.45 


105157.5 


SB(rs)dm;HII 


9.58 


5 


35 


369 


Virgo N Cl. 







12 





TABLE 1 

CONTINUE 



HRS CGCG VCCUGCNGCIC R.A(2000) doc typo K s tot £>(25) vol momb groupAB Other name 













121318 




1-15-10 1 




SAB (s)ab' HII * Sy 










Virgo N CI 








92 


69101 


131 


7255 


3061 




11 


1-10 1-1-1 




S B c ? sp 


10 64 




60 


2317 


Vircro NT CI 






16 


93 


187029 




7256 


4203 


121505 


06 






SAB 0-:; LINER; Sy3 


7 41 




39 


1091 


Coma I CI. 






5 


94 




145 


7260 


4206 


121516 


81 


130126 




SA(s)bc: 


9 39 




10 


702 


Virgo N CI. 








95 


69107 


152 


7268 


4207 


121530 


50 


V ' i' - > .) U.) 






9 44 


l 


96 


592 


Virgo N CI 






07 


96 


69110 


157 


7275 


4212 " 


121539 


36 


i 9 5/1 n^: 




SAc-HII 


8 38 




60 


83 










97 


69112 


167 


7284 


4216 " 


121554 


44 






SAB (s)b:; HII /LINER 


6 52 




12 


140 


Virgo N CI. 






1 1 


98 


69119 


187 


7291 


4222 - 


121622 


52 


131825 


5 


So 


10.33 


3 


52 


226 


Virgo N CI. 




o 


14 


99 




213 


7305 


3094 




00 


13373 1 


.0 




11.25 


° 


93 
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M 58 
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- 


123748.40 
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8 


77 


2 


16 
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Virgo S CI. 
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10 


46 


1 


96 
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Sdm 


11 


76 


1 


89 
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12 




224 
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- 
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33 
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42 
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10 
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96 
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5 
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123922.26 


-053953.3 
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11 
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43 
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V 


rgo Out. 
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11 


LCRSB 123647. 4-052325 
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70204 
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- 
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123948.02 
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11 


11 


1 


10 
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17 
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- 
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151752.1 
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10 


03 


2 


16 
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V 
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1813 
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- 
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7 


46 


4 


76 
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V 
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10 
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- 
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115443.6 
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9 


17 


5 


10 


1645 


V 
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- 
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9 


58 


3 


95 
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1869 


7842 
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- 
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8 


16 


4 


30 
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8 
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2 
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75 


7 


67 
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11 


84 


1 


38 
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9 
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2 


92 
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240 
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21 
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01 
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1939 


7878 
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6 


42 


9 


63 
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11 




251 


15015 




7926 


4666 


- 


124508.59 


-002742.8 


SABc:;HII;LINER 


7 


06 


4 
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1 
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2 
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20 
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12 




263 


129027 




7989 


4725 




125026.61 
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17 


9 


66 
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11 


61 


1 
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77 
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6 
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3 





04 




268 


71060 




8007 


4746 




125155.37 


120458.9 
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2 
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5 
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TABLE 1 

CONTINUE 



HRS 


CGCG 


VCCUGCNGCIC 


RA(2000) 
h m s 


dec 

i » 


type 


K S tot 
mag 


D(25) 


vcl 
km s" 1 


memb 


groi 


ipAB 
mag 


Other name 


271 


100015 
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125244.04 


155055.9 
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10.93 


3 


00 
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0, 
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273 


15031 


8020 


4771 - 


125321.27 


011609.0 
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2 
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12 
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14 
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9.87 


2 


10 
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0. 
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4791 - 


125443.97 
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1 


20 
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2 
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3 


10 
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4799 - 
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S? 


9.89 


1 


60 
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0. 
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8045 
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075434.0 
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91 
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0. 
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4803 - 


125533.67 


081425.8 


Comp 
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50 


2664 
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0. 
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43071 


8054 


4808 - 


125548.94 


041814.7 
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9.04 


2 


.60 
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Virgo Out. 
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3908 


125640.62 


-073346.1 


SB(s)d?;HII 


9.10 


1 


82 


1296 


Virgo Out. 




0. 


15 
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15049 


8078 


4845 - 


125801.19 


013433.0 


SA(s)ab sp;HII 


7.79 


5 


20 


1097 


Virgo Out. 




0. 
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286 


71092 


8102 


4866 - 


125927.14 


141015.8 
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7.92 


6 


00 
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Virgo Out. 




0. 


12 
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8121 


4904 - 


130058.67 


-000138.8 


SB(s)cd;Sbrst 


9.50 


2 


40 
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288 






4941 - 


130413.14 


-053305.8 
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8.22 


3 


63 


1114 


Virgo Out. 




0. 
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289 






4981 - 


130848.74 


-064639.1 


SAB (r) be; LINER 


8.49 


2 


75 


1678 


Virgo Out. 




0. 


18 




290 


189037 


8271 


5014 - 


131131.16 


361654.9 


Sa? sp 


10.11 


1 


70 


1136 


Canes Ven. Spur 


3 


0. 


03 




291 


217031 


8388 


5103 - 


132030.08 


430502.3 


Sab 


9.49 


1 


45 


1297 


Canes Ven. Spur 




0. 


08 




292 


218010 


8439 


5145 - 


132513.92 


431602.2 
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9.33 


2 


00 


1225 


Canes Ven. Spur 




0. 


05 




293 


16069 


8443 


5147 - 


132619.71 


020602.7 


SB(s)dm 


9.73 




91 


1092 


Virgo Out. 




0. 






294 


246017 


8593 
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133601.22 


495739.0 


Sb 


10.42 


2 


19 


1608 


Canes Ven. Spur 




0. 


05 




295 


73054 


8616 


5248 - 


133732.07 


085306.2 


(R)SB(rs)bc;Sy2;HII 


7.25 


1 


79 


1152 


Virgo-Libra Cl. 


P 


0. 


11 




296 


190041 


8675 


5273 - 


134208.34 


353915.2 


SA(s)0;Syl.9 


8.67 


2 


75 


1064 


Canes Ven. Spur 




0. 


04 


KPG391A 


297 


246023 


8711 


5301 - 


134624.61 


460626.7 


SA(s)bc: sp 


9.11 


4 


17 


1508 


Canes Ven. Spur 




0. 


07 




298 


218047 


8725 


5303 - 


134744.97 


381816.4 


Pec 


10.23 





91 


1419 


Canes Ven. Spur 


CP 


0. 


06 


KPG397A 


299 


45108 


8727 


5300 - 


134816.04 


035703.1 


SAB(r)c 


9.50 


3 


89 


1171 


Virgo-Libra Cl. 




0. 


10 




300 


218058 


8756 




135035.89 


423229.5 


Sab 


10.34 


1 


70 


1354 


Canes Ven. Spur 




0. 


06 




301 


17088 


8790 
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135254.46 


-010652.7 


SB(rs)c: 


9.94 


4 


17 


1380 


Virgo-Libra CI. 




0. 


20 




302 


45137 


8821 


5348 - 


135411.27 


051338.8 


SBbc: sp 


10.87 


3 


55 


1443 


Virgo-Libra Cl. 


5 


0. 


12 




303 


295024 


8843 


5372 - 


135446.01 


583959.4 


S? 


10.65 





65 


1717 


Canes Ven-Came. Cl. 




0. 


04 
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46001 


8831 


5356 - 


135458.46 


052001.4 


SABbc: sp;HII 


9.64 


3 


09 


1370 


Virgo-Libra Cl. 
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2 


19 
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13 




306 
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0. 
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46009 


8853 


5364 - 
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SA(rs)bc pcc;HII 


7.80 


6 


76 
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Virgo-Libra Cl. 
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0. 


12 
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46011 


8857 




135626.61 


042348.0 
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11.93 





91 
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Virgo-Libra Cl. 




0. 
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272031 
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11.95 


1 


86 
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Canes Ven-Came. Cl. 
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0. 
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47010 
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035928.4 
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9.98 


3 


72 


1718 
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0. 


13 
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47012 


9175 


5566 - 
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035600.9 
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7.39 


6 


61 


1492 


Virgo-Libra Cl. 


3 


0. 


13 
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312 


47020 


9183 


5576 - 


142103.68 


031615.6 


E3 


7.83 


3 


55 


1482 


Virgo-Libra Cl. 


3 


0. 


14 
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47022 


9187 


5577 - 


142113.11 


032608.8 


SA(rs)bc: 


9.75 


3 


39 


1490 
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3 


0. 
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19012 


9215 




142327.12 


014334.7 


SB(s)d 


10.54 


2 


19 


1389 


Virgo-Libra Cl. 




0. 


14 
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9242 




142521.02 
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11.73 


5 


01 


1440 
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47063 


9308 
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031400.2 
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8.25 


2 


69 
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4 


0. 


14 




317 


47066 


9311 
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143001.85 


034622.3 


S? 


11.70 


1 


10 


1716 


Virgo-Libra Cl. 




0. 


15 
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47070 


9328 


5645 - 


143039.35 


071630.3 


SB(s)d 


9.69 


2 


40 


1370 


Virgo-Libra Cl. 




0. 


12 
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9353 


5669 - 


143243.88 
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SAB(rs)cd 


10.35 


3 


98 
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Virgo-Libra Cl. 




0. 


12 
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47090 
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5668 - 


143324.34 


042701.6 


SA(s)d 


11.71 


3 
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0. 
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S? 
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89 
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S 


9.55 


2 


24 


1640 


Virgo-Libra Cl. 


CP 


0. 
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TABLE 2 
The SPIRE performances 



A 

Ang. Res 
Lin. Res 



250 Aim 
18.1" 
1.7 kpc 



350 fj,m 

25.2" 
2.4 kpc 



500 (im 
36.9" 
3.6 kpc 



Sensitivity 12 mjy/beam 8 mjy/beam 12 mjy/beam 

Note: 1-iT sensitivity (instrumental noise) for a standard map (one 
repeat = two cross-linked scans) as determined during the science 
verification phase (on board observations) . The linear resolution is 
at 20 Mpc. 



4. SPIRE OBSERVATIONS 

With its large field of view (4'x8'), its sensitivity (S ~ 
7 mjy, 5 cr in 1 hour, point source mode), and angular 
resolution (~ 30 arcsec; see Table [2]), SPIRE on Herschel 
is the ideal instrument for the proposed survey. 

In constructing our observing program, we have 
used different integration times for early- and late-type 
galaxies because the former are known to contain less 
dust. 

Integration times for early-types were determined by 
assuming a lower limit of ~ 10 4 M© to their total dust 
mass as determined from IRAS observations (Bregman 
et al. 1998) or measured from optical absorption line 
measurements (Goudfrooij et al. 1994; Van Dokkum 
& Franx 1995; Ferrari et al. 1999). At 20 Mpc, a 
galaxy with a dust mass of ~ 10 4 M© would have a flux 
density at 250 fim of 11 mjy. We estimate that with the 
adopted integration times we should detect an elliptical 
with ~ 10 4 M0 of dust at the 3cr level. 

Integration times for spirals have been chosen so 
that we should be able to detect dust outside the 
optical radiu0, for which there is evidence from ISO 
observations (Alton et al. 1998; Davies et al. 1999). By 
combining ISOPHOT (Alton et al. 1998) and SCUBA 
(Vlahakis et al. 2005) observations of extended sources 
with the spectral energy distribution of normal galaxies 
of different type (Boselli et al. 2003), and assuming a 
standard gas to dust ratio of ~ 100 with the extragalac- 
tic calibration for the dust-emissivity coefficient (James 
et al. 2002), we estimate that the dust associated with 
the extended HI disc would have a flux density of ~ 14 
mjy per beam at 250 fim for HI column densities of ~ 
10 20 atoms cm" 2 . In normal galaxies such gas column 
densities are generally reached well outside the stellar 
disc, at ~ 1.5 optical radii (Cayatte et al. 1994; Broelis 
& Van Woerden 1994). Because of metallicity gradients, 
however, the gas to dust ratio is expected to increase 
in the outer discs, as observed in M31 by Cuillandre et 
al. (2001). We would therefore expect to observe lower 
far-IR flux densities. With 12 minutes of integration 
time (3 scan-maps) we will get to an instrumental noise 
of 6.9 mjy/beam, when the expected emission is 14 mjy 
per beam; this sensitivity will allow us to detect the dust 
emission in the outer disc by integrating flux densities 
along elliptical annuli with a sensitivity of ~ 30 better 
than that of previous observations (Alton et al. 1998). 



All galaxies outside the Virgo cluster will be observed 
in scan-map mode (30" /sec). For each galaxy, the size 
of the scan-map has been chosen based on the optical 
size of the galaxy. For late-type galaxies, the observing 
mode has been chosen so that our images will cover 
the galaxy at least out to 1.5 x D25 (where D25 is the 
25 mag arcsec -2 isophotal diameter) i.e. out to the 
approximate HI radius. For spirals we will make three 
pairs of cross-linked observations, which will reach a 
1 a sensitivity in instrumental noise of 6.9, 4.6 and 
6.9 mjy/beam at 250, 350 and 500 /im respectively as 
determined during the science verification phase. These 
values are comparable to the noise from the confusion 
of faint unrelated sources (4, 5 and 6 mjy/beam in 
the three bands). Elliptical galaxies are likely to be 
weaker and thus need longer integration times, we have 
compromised by only requiring the image to contain 
the galaxy out to D 2 5- For ellipticals and SOs we will 
make eight pairs of cross-linked observations reaching a 
la sensitivity in instrumental noise of 4.2, 2.8 and 4.2 
mjy/beam at 250, 350 and 500 /im respectively, which 
in the latter two bands is below the noise expected from 
the confusion of faint sources. 

Sixty square degrees centered on the Virgo cluster re- 
gion will be covered by Herschel in parallel mode (i.e. us- 
ing both SPIRE and PACS) as part of the Herschel Virgo 
Cluster Survey HeViCS Eighty-three HRS galaxies 
fall into this region. To avoid duplication, these galaxies 
will be observed only during the HeViCS survey. HeViCS 
will survey this region in fast-scan mode (60" /sec) using 
both SPIRE and PACS. The survey will make eight scans 
reaching a 1 cr of instrumental noise of 4.5, 6.2 and 5.3 
mJy/bcam at 250, 350 and 500 respectively, as de- 
termined using the Herschel Observation Planning Tool 
(HSpot v4.2.0). This is not as deep as our observations of 
ellipticals and SOs outside Virgo, but the effect of confu- 
sion noise means that the decrease in effective sensitivity 
is really small. On the plus side, for the HRS galaxies 
in Virgo we will also have PACS images at 110 and 160 
fim. 

5. COMPLEMENTARY DATA 

A number of key aims of our survey (see Section 2) 
can only be achieved through the combination of Her- 
schel observations with corollary data. UV to near-IR 
imaging is needed to determine the distribution and 
content of the different stellar populations, to quantify 
the intensity of the interstellar radiation field and to 
study the recent activity of star formation (UV and 
Ha). Optical spectroscopy is crucial to measure stellar 
and gas metallicities and to determine the presence of 
any nuclear activity. At the same time, the Balmer 
emission line can be used to measure the dust extinction 
within HII regions and are thus essential for quantifying 
the current level of star formation activity. HI and 
CO line data are necessary for determining the content 
and distribution of the gaseous component of the ISM, 
the principal feeder of star formation in galaxies. High 
resolution HI imaging can also be used to study the 
kinematical properties of the target galaxies. Mid-, 
far-IR and sub-mm data, combined with SPIRE data, 
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will be used to study the physical properties of the 
different dust components (PAHs, very small grains and 
big grains), and radio continuum data for measuring the 
thermal and synchrotron emission. 

Given its definition, the HRS is easily accessible for 
ground based and space facilities: the selected galaxies 
are in fact relatively bright (tub < 15 mag) and extended 
(~ 2-3 arcmin). Listed below are the most important 
references for ancillary data. 

5.1. UV, optical, near-IR and Ha imaging 

Of the 323 galaxies in the HRS, 280 have been 
observed by the Galaxy Evolution Explorer (GALEX) 
in the two UV bands FUV (A cff = 1539A, AA = 442A) 
and NUV (A cff = 2316A,AA = 1069A). Observations 
have been taken as part of the Nearby Galaxy Survey 
(NGS; Gil de Paz et al. 2007), the Virgo cluster survey 
(Boselli et al. 2005), the All Imaging Survey (AIS) 
or as pointed observations in open time proposals. A 
GALEX legacy survey has been recently accepted to 
complete the UV observation of the HRS galaxies at a 
deepness of the Medium Imaging Survey (1500 sec/field). 

SDSS (DR7 release, Adelman-McCarthy et al. 2008) 
images in the ugriz filters are available for 313 objects. 
All galaxies have been observed in the JHK filters by 
2MASS (Jarrctt et al. 2003). Deep B,V, H and K 
band images for all Virgo cluster galaxies and for some 
Coma I Cloud galaxies are available on the GOLDMinc 
database (http://goldmine.mib.infn.it/; Gavazzi et al. 
2003). These have been taken with pointed observations 
during the near infrared and optical extensive surveys 
of the Virgo cluster carried out by Boselli et al. (1997, 
2000; 2003) and Gavazzi et al. (2000). 

An Ha+[NII] imaging survey of the star forming HRS 
galaxies found outside the Virgo cluster under way at 
the 2.1m San Pedro Martir telescope is almost complete. 
Combined with images taken in the Virgo cluster region 
(Boselli & Gavazzi 2002, Boselli et al. 2002 and Gavazzi 
et al. 1998; 2002; 2006) Ha+[NII] data are now available 
for 221 of the 258 late- type objects and 26 of the 65 early- 
types. 

5.2. Integrated spectroscopy 
A low resolution (R ~ 1000), integrated spectroscopic 
survey in the wavelength range 3500-7200 A is under- 
way at the 1.93m OHP telescope. In order to sample 
the spectral properties of the whole galaxy, and not just 
those of the nuclear regions (these last provided by the 
SDSS for 106 galaxies in the DR6), observations have 
been executed using the drifting technique described in 
Kennicutt (1992). Exposures are taken while constantly 
and uniformly drifting the spectrograph slit over the 
full extent of the galaxy. A resolution R ~ 1000 
is mandatory for resolving [Nil] from Ha and mea- 
suring the stellar underling Balmer absorption under H/3. 

Data for 64 Virgo cluster galaxies in the HRS have al- 
ready been published in Gavazzi et al. (2004) along with 
a few other galaxies in Moustakas & Kennicutt (2006), 
Jansen et al. (2000) and Kennicutt (1992b). We have 



data from the literature or from our own observations for 
256/258 of the late-types and 33/65 of the early-types, 
making the whole sample complete at 89 %. The re- 
maining objects will be included in the future observing 
runs. 

5.3. HI and CO lines 

Single beam HI observations are available from a 
large variety of sources. Most of the galaxies have 
HI data in Springob et al. (2005) and Gavazzi et al. 
(2005), this last reference limited to the Virgo cluster 
region. Out of the 8 late-type galaxies, 249 have an 
HI measurements, as do 55/65 of the early-types. All 
galaxies in the 0° < dec < 30° range (80 %) will be 
observed during the ALFALFA survey under way at the 
305m Arecibo radio telescope (Giovanelli et al. 2005). 
Data will be gathered with an homogeneous sensitivity 
(rms=2.5 mJy) and spectral resolution (5.5 km s _1 ): at 
the average distance of the HRS the ALFALFA survey 
will detect all sources with MHI > 10 7 5 M©. VLA and 
WSRT HI maps are available for 236 HRS galaxies, from 
the WHISP (van der Hulst 2002) and VIVA (Chung et 
al. 2009) survey, this last limited to the Virgo region, 
from VLA archives or from our own WSRT observations. 

A 12 CO(1-0) line (2.6 mm, 115 GHz) survey of the 
HRS galaxies is under way at the 12m Kitt Peak tele- 
scope. Fifty-eight spiral galaxies have been observed to 
date, with an average rms of 3 mK in T* R scale. Thanks 
to these new observations and to data in the literature 
CO measurements are now available for 161/258 late- 
type and 22/65 early- type galaxies with a detection rate 
of 83% and 55% respectively. Data have been taken 
from different sources, mostly from the FCRAO survey 
of Young et al. (1995), or Kenney & Young (1988) and 
Boselli et al. (1995; 2002) for the Virgo cluster region. 
We obtained time at the James Clerk Maxwell Telescope 
(JCMT) on the 345 GHz HARP array receiver to search 
for CO (3-2) emission from the central 3' x 3' region (with 
15 arcsec resolution) of a subsample of late type galaxies 
with Kstot < 8.7. As of the beginning of January 2010, 
42 of a total of 56 galaxies in the subsample have been 
observed to completion. 

5.4. Mid-IR, far-IR and sub-mm 

Infrared data for the HRS galaxies are already avail- 
able from different sources. 208/258 late-type and 32/65 
of the early-type galaxies have been detected by IRAS 
at 60 and 100 fim, only 103 (Sa-Sd-Im-BCD) and 17 
(E-SO-SOa) at 12 and 25 /im. Virgo galaxies have been 
observed with CAM (45), PHOT (26) and LWS (21) on 
ISO (Leech et al. 1999; Malhotra et al. 2001; Roussel 
et al. 2001; Boselli et al. 2002, 2003; Tuffs et al. 2002). 
Spitzer observations have been completed for 157 HRS 
galaxies with IRAC and 181 galaxies with MIPS. The 
data for all of these galaxies will eventually be available 
from the Spitzer archives. The pipeline-processed IRAC 
data is sufficient for science analysis, but MIPS data has 
been reprocessed using the MIPS Data Analysis Tools 
(Gordon et al. 2005) along with additional processing 
software. 

All HRS galaxies have been observed during the re- 
cently completed AKARI all sky survey in the 9, 20, 
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70, 90 and 160 [im bands which covered ~ 95% of the 
whole sky with an angular resolution comparable to that 
of SPIRE (~ 30 arcsec at 90 fim) and a sensitivity ~ 
three times better than IRAS. Pointed observations with 
AKARI and SCUBA are also available of a few ellipticals. 

5.5. Radio continuum data 

The HRS galaxies have been observed by the NVSS 
survey at 20 cm (1415 MHz) (Condon et al. 1998) with 
an angular resolution of 45 arcsec and a sensitivity of 
2.5 mJy. 22/65 and 159/258 of the early- and late-type 
galaxies respectively have been detected at 20 cm. The 
data of the FIRST survey, at 5 arcsec resolution, are also 
available (Becker et al. 1995). For a fraction of galaxies 
data at 6.3 and 2.8 cm are also available from Niklas et 
al. (1995) and Vollmer et al. (2004). 

5.6. X-ray data 

X-ray data at 0.2-4 keV from the Einstein Observatory 
imaging instruments (IPC and HRI) are available for 31 
early-type and 52 late-type galaxies from Fabbiano et 
al. (1992) and Shapley et al. (2001). Ten HRS galaxies 
have been detected by the ROSAT all sky survey (Voges 
et al. 1999), while Chandra observations of many of the 
early-type objects within the HRS are still under way 
(Gallo et al. 2008). 

Table [5] shows the availability of ancillary data for the 
HRS galaxies. 

6. HRS DATA PROCESSING, DATA PRODUCTS AND WEB 

SITE 

The HRS galaxies are distributed widely across the 
sky, the survey will progress along with the Herschel 
mission, that is expected to last for three and a half 
years, with the first six months for commissioning, per- 
formance verification, science demonstration and the last 
three years for science. Data for each single observation 
will be proprietary for one year after it is taken. Dur- 
ing this time, the data will be processed using the most 
up-to-date pipeline developed by the SPIRE-Instrument 
Control Centre since some members are also part of our 
consortium. This ensures we have access to the latest 
calibration files and bug fixes. Aside from the default 
data processing, the pipeline will also include new steps 
developed by the key-project teams (such as cosmic ray 
removal and baseline subtraction) for producing the most 
accurate representation of large-scale structures, such as 
the extended emission around galaxies. Additional post- 
processing steps will be applied to prepare the data so 
that it can be easily used by the broader astronomical 
community. This will include reformatting the headers 
of the fits files so that they are easy to understand. Fur- 
thermore, the data will have been visually inspected so 
as to ensure that they suffer from no severe artefacts re- 
lated to the observations or data reduction, and in cases 
where problems are identified, the data will be manually 
re-processed to obtain the best results. 
After the proprietary time the team plans to 
make the reduced data available to the community 
through a dedicated web page. Herschel and an- 
cillary data will be archived in an Information Sys- 
tem at the Laboratoire d'Astrophysique de Marseille 
developped under the SITools middleware interface 



TABLE 3 

The completeness in the data sample 



A 
K 
HI 

FIR (60 Atm) 
Radio cont. (20 cm) 



E+SO+SOa 

100% 
85% (40%) 
95% (62%) 
100% (34%) 



Sa-Sm-Im-BCD 
100% 
97% (93%) 
89% (83%) 
100% (62%) 



Note: Values in parenthesis give the detection rate. 

dhttp: / / vds.cnes.fr /sitools / ) with a Postgresql database. 
Images will be stored as fits files, and distributed in fits 
format while catalogues will be distributed in VO table 
and ASCII (cvs) format. The SPIRE data delivered by 
the team will be optimized for use by the wider astro- 
nomical community. 

Given the legacy nature of this project, we will also 
present through an interactive database all the ancillary 
data, X, UV, optical, near- and far-IR, sub-mm and ra- 
dio continuum fluxes and images, optical and line spec- 
troscopy and derived parameters (structural parameters, 
metallicities, SFR...) with the aim of providing the com- 
munity with a unique dataset, a suitable reference for 
future studies. 

7. THE STATISTICAL PROPERTIES OF THE HERSCHEL 
REFERENCE SAMPLE 

In this section, we investigate how representative 
the HRS is of the galaxy population as a whole, thus 
unbiased versus cosmic variance. This is necessary since, 
as selected, the HRS is biased versus galaxies located 
in a high density environment (Virgo). One way to 
test this is to see whether the luminosity distributions 
at different wavelengths in the HRS are similar to the 
local galaxy luminosity functions at the corresponding 
frequencies. 

Shown in Fig. [3] is a comparison of the HRS lumi- 
nosity distribution in the K-band (histogram) to the K 
band luminosity function of Cole et al. (2001) (dotted- 
dashed line; a = 0.96 ± 0.05, M* K = -23.44 ± 0.03) and 
Kochanek et al. (2001) (solid line; a = 1.09 ± 0.06, 
M K = -23.39 ± 0.05), this last also considering sepa- 
rately early- (red dotted line; a = 0.92 ± 0.10, M* K = 
-23.53 ± 0.06) and late-type (blue dashed line; a = 0.87 
± 0.09, M* K = -22.98 ± 0.06) galaxies. The Kochanek 
et al. (2001) values of the K band luminosity function 
have been corrected to take into account the different 
assumed Ho, while the Cole et al. (2001) MJ Kron mag- 
nitudes have been converted into total magnitudes as 
in this work. The zero point of the plotted luminosity 
functions is determined assuming the same number of 
objects as HRS within the given magnitude limit. Fig- 
ure [3] shows that, when considered from the perspective 
of the K-band, the HRS is a good approximation to a 
volume-limited sample down to an absolute magnitude of 
Mk — —20. This absolute magnitude limit is relatively 
high and explains why the HRS is under-represented in 
Im and BCD galaxies (Fig. [SJ. The sample includes only 
the brightest ellipticals and lenticulars, whose K band lu- 
minosity function ends at Mk ~ -18, while it does not 
include any quiescent dwarf system (Mk > -20) . 

Figure [4] shows the HI mass distribution of the HI 
detected HRS early-type (black histogram) and late- 
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type (gray histogram) galaxies compared to the HIPASS 
HI mass function determined by Zwaan et al. (2005) 
(dashed line, a = -1.37 ± 0.03, M H i* = 10 9 - 8±003 
M© and <fi * arbitrarily normalized to match the data). 
The sample includes galaxies with HI masses ranging in 
between ~ 10 s and 10 10 M© but it is clearly under- 
representing galaxies with HI masses in the range 10 75 
M© < M H i < 10 9 M© which are easily detectable by 
an HI blind survey such as ALFALFA. The HIPASS HI 
mass function is dominated by gas rich, star forming 
galaxies, while the HRS also includes quiescent, gas-poor 
early-types and Hi-deficient cluster galaxies, as clearly 
observed in the Virgo cluster (Gavazzi et al. 2005; 2008) 
or in A1367 (Cortese et al. 2008b). The difference be- 
tween the HIPASS mass function and the HRS HI mass 
distribution is due to i) the fact that the HRS is K band 
selected and thus does not include gas rich, low lumi- 
nosity star forming systems with HI masses in the range 
10 7 - 5 < MHI < 10 8 - 5 M© with K band magnitudes 12 < 
K < 16, ii) the presence of cluster Hl-dcficicnt late-type 
galaxies. Indeed the HI mass distribution of galaxies with 
a normal HI gas content (Hl-deficiency H < 0.4) is more 
similar to the HIPASS mass function than the distribu- 
tion of gas poor objects (Hl-deficiency > 0.4). Most of 
the late-type Hl-dcficient objects (Hl-deficiency > 0.4) 
in the HRS are indeed located in the core of the Virgo 
cluster as shown in Fig. [5] Galaxies in the outskirt of 
the Virgo cluster or in the surrounding clouds have a nor- 
mal HI gas content (see Table [4| and can be generally 
considered as unperturbed objects (unless belonging to 
groups or pairs) in the study of the effects of the environ- 
ment on the dust properties of galaxies. Indeed, despite 
the relatively poor statistics, the average Hl-deficiency of 
late-type galaxies in the Coma I Cloud, Leo Cloud, Ursa 
Major Cloud and Southern Spur, Crater Cloud, Canes 
Venatici Spur, Canes Venatici - Camelopardalis Cloud 
and Virgo - Libra Cloud is consistent with that of un- 
perturbed galaxies in the reference sample of Haynes & 
Giovanelli (1984; Hl-deficiency < 0.3). The Hl-deficiency 
of late-type galaxies in the different substructures of the 
Virgo cluster are generally consistent with those observed 
in a larger sample by Gavazzi et al. (1999a). The only 
exception is the Virgo E Cloud that we found dominated 
by Hi-deficient objects. We notice however that the ma- 
jority of the most Hl-dehcient galaxies (Hl-deficiency > 
1.0) of the HRS belonging to this substructure are clas- 
sified as early-types in the VCC. 

The Coma I Cloud, previously thought to include gas 
deficient objects (Garcia-Barreto et al. 1994; Gerin & 
Casoli 1994), is composed of galaxies with a normal gas 
content (Boselli & Gavazzi 2009). A detailed and com- 
plete study of the HI properties of late-type galaxies be- 
longing to the other Clouds has not been done so far. A 
rather normal HI gas content of the spiral galaxies in the 
Ursa Major Cluster, a substructure of the Ursa Major 
Cloud, has been observed by Verheijen & Sancisi (2001). 

30 The Hl-deficiency parameter is defined as as the logarith- 
mic difference between the average HI mass of a reference sam- 
ple of isolated galaxies of similar type and linear dimension and 
the HI mass actually observed in individual objects: HI — def 
= LogM H I re f - LogM H I a i, s . According to Haynes & Giovanelli 
(1984), LogM H I re f = a + b X Log(diam), where a and b are weak 
functions of the Hubble type and diam is the linear diameter of the 
galaxy (see Gavazzi et al. 2005). 



Figure [5] shows the far-IR luminosity distribution of 
the HRS early-type (black histogram) and late-type 
(gray histogram) galaxies detected by IRAS at 60 
fim compared the IRAS 60 /im luminosity function 
determined by Takeuchi et al. (2003) (a = 1.23 ± 
0.04, L* = 8.85 x 10 8 L© and 4> * in arbitrary units). 
The luminosity distribution of the HRS follows the 
luminosity function over two orders of magnitude. The 
low-luminosity cutoff is due to the detection limit of 
the IRAS survey. The lack of very luminous galaxies 
is because Luminous (LIRGs, L 60Mm = 10 11 L©) and 
Ultra Luminous (ULIRGs, Leo^m = 10 12 L©) Infrared 
Galaxies are quite rare in the nearby universe and in 
particular inside rich clusters (Boselli & Gavazzi 2006). 

Figure [7] shows the radio continuum (1.4 GHz) lumi- 
nosity distribution of the radio detected HRS early-type 
(black histogram) and late-type (gray histogram) galax- 
ies, as compared to the 2dF/NVSS luminosity function 
of Mauch & Sadler (2007) (solid line). This radio lumi- 
nosity function includes the contribution from both AGN 
(dashed line) and quiescent (dotted line) galaxies. The 
radio luminosity distribution of the HRS is quite typi- 
cal of normal, nearby galaxies (10 20 < Ljiadio < 10 22 W 
Hz" 1 ; Condon et al. 2002). The two brightest sources 
are the radio galaxies M87 (Virgo A, Ljiadio — 10 24 8 
W Hz" 1 ) and M84 (LR adlo = 10 23 3 W Hz" 1 ). The ra- 
dio continuum luminosity distribution as a whole agrees 
well with the field luminosity function, indicating that 
the contamination of cluster galaxies with an enhanced 
radio continuum emission is minimal (Gavazzi & Boselli 
1999c,d). 

In summary, since the K-band light traces stellar mass, 
the good agreement between the HRS luminosity distri- 
bution and K-band luminosity function implies that the 
HRS can be treated as a volume-limited sample down to 
a minimum stellar mass. Figure [5] shows that the HRS 
is also representative of a far-IR selected sample once 
the most active and rare in the local universe LIRG and 
ULIRG are excluded. The HRS matches also the dis- 
tribution of a radio continuum selected sample, although 
the statistic for the brightest radio galaxies is poor. How- 
ever, Fig. UJshows that when viewed from the perspective 
of gas mass, the HRS is not a fair sample of the local uni- 
verse. The presence of Hi-deficient objects makes it ideal 
for studying the effects of the cluster environment on the 
dust properties of nearby galaxies. 
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Database (NED) which is operated by the Jet Propul- 
sion Laboratory, California Institute of Technology, un- 
der contract with the National Aeronautics and Space 
Administration, and of the GOLD Mine database. A.B. 
wishes to thank S. Boissier for his help in the construc- 
tion of the HRS database. We are grateful to the anony- 
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TABLE 4 

The average HI-deficiency of late- type galaxies in the different substructures 



Substructure average standard deviation N. of objects 



Virgo A 


n 01 
u.y 1 


4-0/17 


OO 


T T: T~> 

Virgo B 


0.71 


± 0.46 


25 


Virgo N Cloud 


0.39 


± 0.24 


12 


Virgo E Cloud 


0.83 


± 0.53 


12 


Virgo S Cloud 


0.47 


± 0.47 


26 


Virgo Outskirts 


0.26 


± 0.60 


38 


Coma I Cloud 


0.07 


± 0.35 


8 


Leo Cloud 


0.11 


± 0.31 


45 


Ursa Major Cloud 


0.05 


± 0.34 


15 


Ursa Major Southern Spur 


-0.02 


± 0.20 


4 


Crater Cloud 


0.13 


± 0.40 


4 


Canes Ven. Spur & Camclopardalis 


0.19 


± 0.03 


9 


Virgo - Libra Cloud 


0.31 


± 0.35 
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Note: Upper limits are considered as detections 
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Fig. 1. — The sky distribution of the HRS for early-type (E-SO-SOa; circles) and late-type (Sa-Sd-Im-BCD; crosses) galaxies (left panel). 
Dashed contours delimit the different clouds. The big concentration of galaxies in the centre of the figure is the Virgo cluster (red colour) 
with its outskirts (dark green). Orange symbols are for Coma I Cloud, magenta for Leo Cloud, brown for Ursa Major Southern Spur and 
Cloud, cyan for Crater Cloud, light green for Canes Venatici Spur and Camclopardalis and blue for Virgo-Lybra Cloud galaxies respectively. 
The Virgo cluster region (right panel): black contours show the diffuse X-ray emission of the cluster (from Bohringer et al. 1994). Red 
symbols are for galaxies belonging to the Virgo A cloud, blue to Virgo B, orange to Virgo E, magenta to Virgo S, cyan to Virgo N and 
dark green to the Virgo outskirts, as defined by Gavazzi et al. (1999a). 
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Fig. 2. — The distribution in morphological type of the HRS for 
all (empty histogram) and cluster (filled histogram) galaxies. The 
cluster sample is composed of galaxies members of the Virgo A and 
B clouds. 
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Fig. 3. — The K band luminosity distribution of the HRS for all 
galaxies is compared to the 2MASS K band luminosity function 
of Kochanek et al. (2001) (black solid line) and Cole et al. (2001) 
(black dotted-dashed line) (upper panel). The K band luminosity 
distributions of the E-SO-SOa (central panel) and Sa-Sd-Im-BCD 
(lower panel) galaxies in the HRS are compared to the Kochanek 
et al. (2001) K band luminosity function of early-type (red dotted 
line) and late- type (blue dashed line) galaxies. Poisson errors are 
also indicated. 
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Fig. 4. — The atomic hydrogen mass distribution in linear (upper 
panel) and logarithmic (middle panel) scales for the HI detected 
HRS galaxies (empty histogram). The black and gray histograms 
are for early-type (E-SO-SOa, 35% detected) and late-type (Sa-Sd- 
Im-BCD; 93 % detected) galaxies respectively. Lower panel: The 
atomic hydrogen mass distribution for late-type galaxies with a 
normal HI gas content (Hi-deficiency < 0.4; tiny spaced hashed 
histogram) and gas poor galaxies (Hi-deficiency > 0.4; large spaced 
hashed histogram). The red dashed line is the HI mass function 
determined by Zwaan et al. (2005). The vertical, dashed line 
indicates the detection limit of the ALFALFA survey (2.5 mjy) 
at the distance of 20 Mpc. 
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Fig. 5. — The variation of the Hi-deficiency parameter as a 
function of the angular distance from the core of the Virgo cluster 
(M87). Red circles are for Virgo cluster galaxies (Virgo A, B, N, 
S and E), dark green circles for galaxies in the Virgo outskirts, 
orange triangles for the Coma I Cloud, magenta crosses for the Leo 
Cloud, brown hexagons for Ursa Major Southern Spur and Cloud, 
cyan three skeletal for the Crater Cloud, light green squares for 
Canes Venatici Spur and Camelopardalis and blue stars for the 
Virgo-Lybra Cloud. Empty symbols are for Hi-detected galaxies, 
filled symbols for upper limits. 
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Fig. 6. — The infrared luminosity distribution in linear (upper 
panel) and logarithmic (lower panel) scales for the FIR (60 fim) 
detected HRS galaxies (empty histogram). The black and gray 
histograms are for early-type (E-SO-SOa) and late-type (Sa-Sd- 
Im-BCD) galaxies respectively. The dashed line is the IRAS 60 
fim luminosity function determined by Takeuchi et al. (2003). The 
vertical, dashed line indicates the typical detection limit of IRAS 
(0.4 mjy) at the distance of 20 Mpc. 



Fig. 7. — The radio continuum (1415 MHz) luminosity 
distribution in linear (upper panel) and logarithmic (lower panel) 
scales for the detected HRS galaxies (empty histogram) . The black 
and gray histograms are for early-type (E-SO-SOa) and late-type 
(Sa-Sd-Im-BCD) galaxies respectively compared to the 2dF/NVSS 
luminosity function of Mauch & Sadler (2007) (solid line). This 
radio luminosity function includes both the contribution of AGN 
(red dashed line) and quiescent (blue dotted line) galaxies. The 
vertical, dashed line indicates the typical detection limit of NVSS 
(2.5 mjy) at the distance of 20 Mpc. 
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Ancillary data for the Herschel Reference Sample. 
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TABLE 5 

CONTINUE 



2MASS SDSS 



IRAS 20cm 
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